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Abstract 

We link the recent computation beyond leading twist of the impact factor of 
the transition 7^ — > px performed in the light-cone collinear approach, to the 
dipole picture by expressing the hard part of the process through its Fourier 
transform in coordinate space. We show that in the Wandzura-Wilczek ap- 
proximation the impact factor up to twist 3 factorises in the wave function 
of the photon combined with the distribution amplitudes of the p— meson 
and the colour dipole scattering amplitude with the t— channel gluons. We 
show also that beyond the Wandzura-Wilczek approximation, the hard con- 
tribution of the amplitude still exhibits the signature of the interaction of 
a single colour dipole with the t— channel gluons. This result allows a phe- 
nomenological approach of the helicity amplitudes of the leptoproduction of 
vector meson, by combining our results to a dipole/target scattering ampli- 
tude model. 
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1. Introduction 

In the limit of asymptotical energy, the understanding of hadronic scat- 
tering processes is a longstanding question, which can be addressed in per- 
turbative QCD when a hard scale, generically denoted as Q 2 , justifies the ap- 
plicability of perturbation theory. In this so-called perturbative Regge limit, 
in which s 3> — t ~ Q 2 3> A.q CD , the scattering amplitude is dominated by 
the exchange of reggeized gluons (named gluonic reggeons) in the t— channel. 
This regime is governed by the BFKL equation, derived in the leading order 
approximation (LLA) in Refs. [H, Q, S, 0] • The scattering amplitude is a con- 
volution in transverse momentum of impact factors - describing the coupling 
of the targets with reggeons - with the reggeon-reggeon scattering amplitude, 
in the /cy-factorisation framework. 

In the same kinematics, a dipole representation M was proposed, and its 
dynamics was studied at LLA in Refs. \^ and 0, la]. In this picture, the 
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degrees of freedom are coloured dipoles. The scattering of two colourless ob- 
jects (heavy meson for example, or virtual photon, called generically onium) 
is then described completely in terms of these dipoles, and does not involve 
the notion of reggeization. It was checked at the level of the kernel ^ and 
then at the level of full amplitude for onium-onium scattering 10[ that both 
BFKL and dipole models are equivalent, at large N c . 

At even higher energies, based on the growing of partonic densities which 
should be limited due to unitarity arguments, many theoretical investigations 
have been made in order to resum perturbative sub-series which contribute 
to the scattering amplitude, which could be responsible for the required non- 
linearities leading to recombination effects between partons. This includes 
the Generalized Leading Log Approximation, which takes into account any 
fixed number n of t-channel exchanged reggeons [ill, 12, 13, 14 , and the Ex- 



tended Generalized Leading Log Approximation (EGLLA) [15|, |16|, |17J, lla [19 
in which the number of reggeon in t— channel is not conserved. In EGLLA, 
the simplest new building block is the triple Pomeron vertex 17], [l6|, UM, l20| . 



In the Wilson line formalism, non-linear equations have been derived [2ll . 
based on the concept of factorisation of the scattering amplitude 



22, 23, 24 



in rapidity space and on the extension of the operator product expansion 
technique to high-energy Regge limit. Its simplest version, the Balitsky- 
Kovchegov (BK) equation, has also been derived independently by Kovchegov 



25l [26[ in the dipole model. The triple Pomeron vertex is the building block 
responsible for non-linearities in the BK equation. Its complete expression 
beyond large N c , known in EGLLA 
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line formalism in 27 



was rederived based on the Wilson 
in a very compact way. Simi lar l y t he Wilson line 



approach, the Color Glass Condensate (CGC) |28|, |29|, |30| M, |32|, |33|, [3J, [35 
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leads to an equivalent set of non-linear equations. 
The dipole model itself has been the basis of many studies in order to 
unitarize the theory jsl, 0, [38| . As for other models, it can describe both the 



Pomeron and the Odderon |39j] degrees of freedom [40| . Besides its dynamics 
at large s, the dipole representation of the probe has been used to build ob- 
servables sensitive to saturation effects. In particular, the geometrical scaling 
is a natural consequence of saturation. Although strictly speaking, geometri- 
cal scaling does not implies saturation, the fact that geometrical scaling has 
been seen at HERA for Deep Inelastic Scattering (DIS) |4lj], for moderate Q 2 
(virtuality of the photon probe) and very small £ is a good sign that QCD 
probes the saturation regime. Basically, geometrical scaling means that the 
7*P total cross-section, which is a function of Q 2 and x, can be described as a 



3 



function of a single variable, which is ratio of Q and of a saturation scale Q 
In the dipole representation, a simple model was introduced by Golec-Biernat 



for describing the total 7*p cross-section [42 



and Wiisthoff (GBW) [42|, |43 
as well as diffractive events 43] . The basic idea is to describe the 7* — p 
interaction at small x as the scattering of a qq pair (a dipole), formed long 
before the scattering off the proton (in a non-symmetric frame where the nu- 
cleon is at rest). This initial dipole is characterized by a transverse size r and 
by a relative fraction a of longitudinal momentum carried by the quark and 
the antiquark. One should then parametrise the dipole-nucleon scattering 
cross-section. The total cross-section can be expressed as (here r = \r\) 



°t,l{ x 'Q 2 ) = fd 2 L^da\ty T ,L(a,r)\ 2 &(x,r 2 ) 



(1) 



where ^t,l is the photon wave function for the transverse (T) and longi- 
tudinally polarized (L) photons, which can be expressed in terms of Bessel 
functions. One then uses an effective description of the saturation dynamics, 
implemented in a simple functional form for the effective dipole cross section 
a(x,r) which encodes the interaction of the qq dipole with a nucleon, 



2\ 



a x, r 



j. 2 



where the ^-dependent radius Rq is given by 



1 / x \ A/2 

^o(x) = —[-?-). (3) 



GeV \x J 

One gets for a(x, r 2 ) a usual linear perturbative description in the regime 
r R Q . In that limit the dipole cross-section indeed scales like r 2 , a typical 
behaviour for colour-transparency. This is for example what is obtained when 
computing the dipole-dipole scattering through a two gluon exchange (form 



of which just comes out of an eikonal treatment [10]), in a simple model 
where the nucleon probe is made of dipoles. The power-like behaviour of 
Ro{x) is inspired by the hard Pomeron description of the HERA data in the 
linear regime, obtained within the Mueller dipole model when dressed by 
small-x gluon emissions a la BFKL (dipoles in Mueller approach) the initial 
qq dipole. 
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On the other hand, when r gets larger, of the order of 2R , the cross- 
section saturates, with for small Q 2 a maximal value which equals o"o . The 
critical saturation line is thus given by Q = Q s {x) with 



Q s {x) = l/R {x). 



(4) 



To get the geometric scaling prediction, one should note that the non-trivial 
functional dependency of the 7* wave function is completely expressed as a 
shape in the dimensionless variable f = \r\ Q . In this variable, the dipole 
cross-section (TSl) has now the functional form 



a 



QM 
Q 2 



(5) 



after using the saturation scale 01]). Combining Eqs. (0) and (JTJ, it is now 
clear that the functional form of the dipole cross-section (J2J) remains intact 
after the perturbative dressing due to the 7* coupling, based on the fact that 
the squared wave functions are peaked at r ~ 2/Q . One can therefore expect 
the geometric scaling of the total cross-section at small x 



a 



tot 1 



a 



tot 



Q 2 /Ql 



,x 



(6) 



This saturation model was further extended by including the effect of DGLAP 



evolution 

„2 



44 



This can be done by replacing in the dipole cross-section 
a[x, V) the elementary r 2 dipole-dipole cross-section of the colour trans- 
parency regime by a more elaborate cross-section which includes DGLAP 
evolution through the gluon density entering the scattering off the dipole. In 



Ref. 45 



46, 47, 48, 49, 50 for F 2 has 



a compilation of all available data 
been used, leading to an almost perfect agreement with geometric scaling. 

Another sign of the saturation effect was presented in Ref. 5l| were a 
model, in the spirit of GBW, was constructed inspired by the BK equa- 
tion, smoothly interpolating between the saturated regime and the linear 
regime. This model was able to describe the most recent very precise F 2 
data 46|, |47J, |48| in the moderate domain of Q 2 , explicitly exhibiting a need 



for saturation effects at small x and small Q 2 in comparison with a pure linear 
evolution a la BFKL. Geometrical scaling can be extended at each impact 
parameter b, as shown in Ref. 52] . The idea of a b— dependent saturation 



scale is due to Mueller [53[, while a b— dependent S— matrix for dipole-proton 



scattering was studied for diffractive meson electroproduction in Ref. [54 
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Beside this framework whose main focus is devoted to QCD dynamics 
in the asymptotical limit, there has been much progress in the understand- 
ing of the partonic structure of mesons at leading and beyond leading twist 
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This is mainly based on exclusive vector meson production at inter- 
mediate energies (HERMES, JLab) and asymptotical energies (HI, ZEUS). 
The natural theoretical framework is collinear factorisation, which describes 
the scattering amplitude as a convolution (in the space of longitudinal mo- 
mentum fraction) of an hard part with a Distribution Amplitude (DA). This 
framework is under control at leading twist for many processes, but faces 
severe potential divergencies, in particular for processes involving or domi- 
nated by contributions beyond leading twist. Recent data have shown that 
some 7* — > p transitions, which vanish at twist 2, are copiously produced 



and thus deserve an understanding from first principles. In Refs. |56l . 157 
a new framework was developed in order to deal with such transitions, and 
have been recently applied to describe HI and ZEUS data [58|. In this study, 



a non negligible contribution from the low vituality t— channel exchanged 
gluons was found. This could be a sign that saturation effects are important 
for both twist 2 and twist 3 7* — > p transitions. 

The natural question which we want to address in this paper is there- 
fore whether it is possible to construct a consistent framework to deal with 
exclusive processes, for example exclusive vector meson electroproduction, 
beyond leading twist, and at asymptotical energies where saturation effects 
are expected. Our strategy is to combine the dipole framework with collinear 
factorisation beyond leading twist, that is to construct the 7* — > p transition 
in the coordinate space for the t— channel gluons, for finite N c , including 
both 2-parton and 3-parton contributions. This treatment is different in 
spirit from the pure 2-parton approach based on models for the light-cone 



p— meson wave function [59|, [60|, |6l|, [62|, |63J. Our result shows that the cou 



pling of the proton to this transition only involves a dipole coupling, and 
does not require any quadripole structure, which in principle can arise in 



high energy scattering [38] . The phenomenological implementation of our 
result, which requires to combine it with a model for the proton - dipole 
scattering amplitude, is left for a forthcoming publication. 

The structure of the paper is the following. In Section [2J we derive the 
dipole representation of the twist 3 7* — > p impact factor in the 2-parton 
approximation. The crucial ingredient to get this dipole representation is 
the use of the QCD equation of motion (EOM). In Section [31 we extend this 
analysis to the 3-parton sector, again using the QCD EOM including now 
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Figure 1: a: 7* — > p impact factor, b: k integration contour entering the definition of the 
7* — > p impact factor. 



genuine twist 3 effects. This leads to the main results of our analysis, shown 
in Eqs. ( 1158111591) . 

2. 2-parton impact factors up to twist 3 in the dipole picture 

2.1. Impact factor representation 

In the kx factorisation approach |64|, |65|, |66|, |67j, the description of the 
scattering amplitude for the process 

l*(q) + N(p 2 )^p( Pl ) + N (7) 

involves the 7* — > p impact factor $ of the subprocess 

g(k u £1) + 7*(g) -)• g{k 2 , e 2 ) + p(pi) ■ (8) 

This impact factor is illustrated in Fig. [T^,, in the kinematical region where 
virtualities of the photon, Q 2 , and of t— channel gluons k\, are of the same 
order, Q 2 ~ kj_, and much larger than A.q CD . Neglecting masses, one consid- 
ers that in the reaction (J7J) the vectors p\ and p 2 are light-like. This sets a 
natural Sudakov basis for momenta, and we denote 2p\-p 2 = s. The 7* — > p 
impact factor is defined as the integral of the S"-matrix element S]^ g Pt9 
with respect to the Sudakov component of the t-channel k momentum along 
p 2 , which is directly related to k = (q + k\) 2 , the Mandelstam variable s 
for the subprocess (|S]), as illustrated in Fig. [IJi. Closing the k contour of 
integration below shows that the impact factor can equivalently be written 
(see Fig. [Tb) as the integral of the K-channel discontinuity of the S'-matrix 
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element^ Sf 9 ^ pT 9 



+00 

$ r ^ p (k,r-k) = e 7> — / — S^ 9 ^ p9 (k, r-k) 



-00 
+00 



ol A 4 



dn 

s J 2tt 





Disc K S^ 9 ^ p9 (k } r-k) 



(9) 



The two reggeized gluons have so-called non-sense polarisations E\ 




+ 




Figure 2: 2- and 3-parton correlators attached to the hard scattering amplitude of the 
7* — > p impact factor, where vertical lines are hard t— channel gluons in the colour singlet 
state. 

Up to twist 3, the impact factor receives contributions both from 2-parton 
and 3-parton diagrams, as illustrated in Fig. |5J In this section we only con- 
sider the 2-parton contributions, in the so-called Wandzura-Wilczek (WW) 
approximation. The 3-parton contribution will be discussed in Section [3j 

2.2. Collinear factorisation beyond leading twist, in the transverse coordinate 
space 

The impact factor contains hard and soft parts, illustrated in Fig. [21 that 
we need to factorise out in Dirac, colour and momentum spaces. Using the 



Fierz identity in the standard form [68 



X = x n T c 



^r Q Tr (Xr a ) = Tr (XT C 



(10) 



1 In this paper, underlined letters like I denote euclidean momenta, such that l\ = — . 
Overlined letters like y denote 1 — y. 
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with r r 



(r 



a\-l 



the contribution of the 2-parton Fock state to the ampli- 



tude defining the impact factor 7* — > p reads 



7 ->p 



d: [ 



(27T)' 



#(£) S(£) 



d 4 



4 J (2tt) 4 



ti(H(l)Y a ) S 



:n) 



as illustrated in Fig. [3j Spinor indices, as well as colour indices by using Fierz 
identity in colour space, of the hard and the soft parts are now disconnected, 
which achieves the factorisation of the amplitude in Dirac and colour spaces. 





Figure 3: Factorization with Fierz identity of the 2-parton 7* — > p impact factor. 



Using Sudakov variables for the loop momentum £ = ap + fin + £j_ and 
for its Fourier conjugate coordinates z = a z p + /3 z n + z±, the £ integration 
reads 



1 / dy 



x 



d*ze- u -*{p(p)$(z) r a ^(0)|0) 
x (p(p)|^(An + ^)r a ^0)|0). 



dap f dfit 
2tt / 2tt 



,iA(« f -j/) 



27T 



— ( 

2vr 



(12) 



We will denote for conciseness H ra (y,£±) = tr[H (y, £±)T a ] the hard part of 
the amplitude and we denote H ra (y,x±) its Fourier transform in the trans- 
verse xt plane. The amplitude rewritten in terms of H ra (y,x±) reads 



d 2 x± H ra (y, x±) e 



-ix±-£± 



d\ 
( 

2vr 



-i\y 



X 



d 2 z±e 



-it 1 -z 1 



<p(p)$(An + * ± )r«V(0)|0). 



(13) 



2 The minus sign arises due to the fermionic nature of quark and antiquark fields. 
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—e- iXy ( P (j>m\n) 1 »mn = rn p f p [ Vl {y) (e* p ■ n) p, + <p 3 (y) e^J ,(15) 



The collinear approximation corresponds to a Taylor expansion of the term 
e -i£j_-x ± arounc [ fag collinear direction, i.e around £± = 0. Up to twist 3, we 
need to expand up to the first order: e~ te± ' x± = 1 — i£± ■ x± + o(twist 3) . We 
denote by $2-parton the contribution of the k th term in the Taylor series. 
The first term of the Taylor expansion gives terms of twist 2 and 3, as 

= ~\ JdyJ dWH r ly,x ± )J^e-^(p(p)m\n)T a m\0) • (14) 

The corresponding matrix elements of the relevant non-local correlators with 
a light-like separation between quar k fields are parametrised by a set of dis- 
tribution amplitudes (DAs) j56], 57 

r dX 
2^ 

^e-ttf(p(p) |VS(An) 757/, V>(0)|0> = m p f p i <p A (y) e paM e*«p? . (16) 

The second term of the Taylor expansion gives only terms appearing at twist 3 
and reads 

x j d 2 z ± e- ll ^(p(p)\tP(Xn + z ± )T a 4j(0)\0) 

x J d 2 z ± ^(e- ie ^)(p(p)m\n + z ± )T^(0)\0) ■ (17) 

An integration by parts, denoting d^= \{d^ — d^), and using the identity 
- -^(P(p)\^n + z ± ) I> V(0)]0> = (p(p p )mz)Te ^(0)|0> , (18) 
finally leads to 

1 

4 

^e-^(p(p p )m\n)T^m\0}. (19) 
10 



<2£L = -7 Idyl d 2 x ± xlH r \y,x ± ] 



The corresponding matrix elements of the relevant non-local correlators in- 
volved in Eq. f[T9j) are parametrised by [56|, [57 



^e-^ipipMXn)^ m\0) = ^m p f p ^(y)p,e; Ta , (20) 

^e- a y(p{p)\i;{\n) l5l , ^(0)|0) = m p f p y T A {y) p p e a ^ e^p & n< . (21) 

To summarize our approach, the impact factor up to twist 3, after performing 
the collinear approximation, has the form 

i 



k=0 

X E ( ~ X± k \ d±)k ( f ^- lX y(p(p)\i;(Xn + z ± )T,m\0)) -(22) 
k=o ' ^ ' z±=o± 



Using the hard and soft parts together, we thus obtain, for the vector part 
of the 2-parton impact factor $y ~* p and axial- vector part <§P A ^ p 

$$^ p =~m p f p JdyJd 2 x x [pz{y)e* p±p - i^(y)p lp (e p± ■ x x )]H r {y, x±), (23) 

®A -\ m pfp M [i<PA(y)£ne* . Pl n + ¥A(y)Plp £ x±e* . Pl n]H^ (y, X±) . 



(24) 

In the following part we will compute the 2-parton hard parts H 1 ^ (y , t j_) 
and H 15 ^ (y,£±), and then we will derive the expressions of their Fourier 
transforms in the transverse coordinate space H 1 '*(y,x±) and H 151 ^ (y , x ±) . 

2.3. Calculations of the hard parts in transverse coordinate space 

We now compute the hard parts in momentum space if 7 ^?/, £j_) and 
H 151,i (y,£±) in a kinematics extended to the transverse degrees of freedom 
£± of the partons when compared to the kinematics of Ref. (56), 57, 58]. 
The partons are kept on the mass-shell and in the collinear limit (£± — > 0), 
this kinematics is the same than in the previous references, which allows to 
compare the final results after integration over x of the Fourier transforms. 
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The Sudakov decomposition of the momenta ara^l 



f 

?q = VPl + £q± H P2 

ys 

e 

£q = VPl ~ V + — P2 

ys 

I 2 1 p 2 
Pp = Pi + - — -P2 =Pi + —P2 • (25) 
s yy s 



The momentum of the incoming photon is 



s 

while the momenta of the gluons in t-channel are 



q=Pi- —P2 , (26) 



K + k 2 + Q 2 
h = p 2 + k± , 

s 

k + k 2 — p 2 

k 2 = = ^p 2 + k ± . (27) 



We use the same labelling of diagrams as in Ref. |57j. After computing 
all the 6 diagrams (a), (b), (c), (d), (e), (f), shown in Fig. HI we perform 
the integral over k by the method of residues to get the contribution to the 
impact factor, according to Eq. (Q. Four poles in k appear, which are 



Kl = 


(l-ykf 

yy 


irj 




K 2 = 


(i+yk) 2 

yy 


IT) 




«3 = 




Q 2 - 


i + y(k 2 + Q 2 ))+ir} 


K4 = 


f((k- 


2 - 


I 2 + y(k 2 + Q 2 )) + irj 



The amplitudes associated to each diagrams (a), (b), (c), (e), that will 
give a non-zero contribution when we integrate over k, closing the con- 
tour in the lower k plane, denoted by C_, reads for the different structures 
Y a = {7 M ,7 M 7 5 } 



3 In this paper, overlined letters like y denote 1 — y. 
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(b) (c) (d) 





(e) (f) 

Figure 4: The 6 hard diagrams attached to the 2-parton correlators, which contribute 
to the 7* — > p impact factor, with momentum flux of external line, along pi direction. 
These drawing implicitly assume that the two right-hand side spinor lines are closed on 
the various possible Fierz structures T, here j) or ^7 5 . 



diagram (a) 



iM v ; - lCh 



diagram (b) 
iMF = 



2ieg 2 5 ab tifai-yfa+MfcT"] 

y/2 2N C (f + ^ 2 )(«-«i) ' 1 ' 



le g 2 s ab 2 tr[Myfc + h + + h + h)fcr a ] 



diagram (c) 



y/2 2N C s yy(n - k 2 )(k - k 3 ) 

2ie g 2 5 ab 1 
a/2 2N c yy(n - k 2 )(k - k 3 ) 

x {ytrfaifa + - ytr(Mh + h)h^ a )} ■ (29) 



ra 2ieg 2 6 ab Jr[My]/>i + h)^ a } 
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diagram (e) 

lie g 2 5 ab 



iMi 



y/2 2N C yy(n - ki)(k - k 4 ) 

x {ytr(Mh -hW) -ytr(Mh -k)h? a )} -(31) 



In the above expressions, the scale fi 2 is defined by 

fi 2 = yyQ 2 . (32) 

The impact factor is proportional to the integral over k of the hard part, 
according to the definition (jHJ), which we symbolically write 

h L ^ - -h^Mn 03) 

with iM va = ELi iM T is the 

sum over the 4 contributing diagrams, and 
since the integration is performed along the contour C_, the r.h.s of ( )33|) 
means the sum of the residue at K\ and ki- We note that 

Ki-« 4 = — ((l-k) 2 + fi 2 ), (34) 



2/2/ 

k 2 -k 3 = — ((/f+A;) 2 + ^ 2 ). (35) 
2/2/ 



Thus, 



e,? 2 ^ 6 / 2/tr[^(^ - y/^jfeF*] - y tr^ 2 (^i + h)^ c 



V22N c s I f + /i 2 

{£ - k) 2 + [i 2 
(£ + k) 2 + H 2 

The vector and the axial-vector hard part read respectively 



e^-l e^-d + k) e^H-k) 



2 fidb yy e g2 



f + /i 2 (l + k) 2 + fi 2 (£-k) 2 + fi 2 
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(36) 



and 



.eg 2 5° 



x e 



fll/pcr 



sy/22N c 



+ /1 2 



{£ + k) 2 + fi 2 



{£ - kf + fi 2 



(37) 



The Fourier transforms of propagators in ( 1361 1371) are related to the modified 
Bessel functions K v (x) 



1 



r + 



2tt 



K (/x|x|)e^ 



+/i 2 



(38) 
(39) 



The Fourier transforms of the vector and axial- vector hard part read thus 
eg 2 5 ab 



(y,x) = 4- 



{-yyKMx\)e!t 



+ p%(y - y)iJ^KMx\) [(1 - e**)(l - e"**) - l] 



(40) 



H ri5 (y,x) = 
eg 2 5 ab 
s(2tt)V2 2N c 



e^ pa e 



7 !/ 



\X\ 



P2v 



[(! _ e **-*)( 1 _ e -t*-s) _ x]. ( 41 ) 



TTie 2-parton 7^ — > impact factor 

Substituting in the expressions of the impact factors, Eqs. ( l23l l24l) . re- 
spectively the vector (1401) and axial-vector (1411) hard parts, we get 



C afe Q 2 



f d 2 x 

dy / 2^{-2^3(?/)-^o(^|x|)e T -e* 



+ (V - V)tf(vM ■ *) ^-liKMd) [(1 - - e"^) - l][ (42) 



and 



-C a6 Q 2 



d 2 x± 
2vr 



x ^i(/i|x|) ((1 - e* 2 )(l - e"^) - 1) 



(43) 
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where we define 



C 



ah 



eff 5 ab 1 



(44) 



Note that the term with if a in Eq. (|24|) vanishes due to the structure of the 
expression flU). Using £x ± e* xPl p 2 £xj_e 1 _ LPiP2 = x( x l e 7-L ' e p_L _ ( e 7-L ' ' 
e*j_)), the axial- vector contribution takes the form 



C ab Q Q 



d 2 x 



d y J ■^■ ( PA.(y) 

x/i|x|i<: 1 ( A t|x|)((l-e^)(l-e- 
The whole 2-parton contribution thus reads 



)-!)• 



5,7 ~>P 
2— parton 



iab/ 



(45) 
(46) 



dy {-2yy<p 3 (y)K (iJ,\x\) e y ■ e* 



+ 



((y - y)<fl(y) - va(v)) 



\x\ 2 



x /i|x|A\(/i|x|)((l - e^ £ )(l - e~ 1 ^) - 1)} 



(47) 



This result does not seem to be proportional to the familiar dipole factor 



N{x : k) 



ik-x 



)(1 



-ik-x\ 



(48) 



describing the coupling to the two t— channel gluons. The dipole structure 
( |48|) is recovered by the use of the QCD EOM, which relate the various DAs 
appearing in Eq. (|47|) . The solution for a DA ip{y) of these equations can 



be split in f(y) = ip (y) + ip 9en (y), where </? (y) is solution of EOM in 
the Wandzura-Wilczek approximation, i.e. in the limit of vanishing 3-parton 
DAs. These equations j56|, 5?J lead in the WW approximation to the relation 



2yy V >Y w (y) + (y-y)tf ww 

which appears in Eq. ( 14T1) rewritten in the form 

C ab Q 2 



+ <p T A ww 



(y) = o 



<f)7 

2— parton 



. d 2 x 



(49) 



(50) 



x { [~(2y y cp 3 (y) + {y~ y)<fl + (Pa(v))] K o(vk\) e 7 • e* 



+ 



[(y - y)fi(y) - fliv)] 



-P — — -7 1 TV \ * 

+ ^(?/)e p -e. 



x /i|xJifi(/x|xJ)A/"(x, &)}, 
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in which we used the following integral relations 

^IxlK^x]) = 2 J d -^KMx\) = j 2 (51) 
f^£-x^t*KM*b) = ■ e* p J ^=K (^\x\) . (52) 

2. 5. The Wandzura- Wilczek approximated impact factor in the dipole picture 

Using Eq. (]49p . the 2-parton WW contribution (150]) exhibits the interac- 
tion term (1 — e*--)(l — e~ 1 --) between the dipole of transverse size x and 
the gluons in t-channel of transverse momenta k as 

x ( -^[ 2 -'^ + V T A ww (y)e; ■ e\ ^\x\KMx\) A%, k) . (53) 

This result will be extended in part 14.11 to the genuine DAs ip gen solutions 
of EOM in the presence of the 3-parton DAs. Similarly to the momentum 
space analysis, one can split the result f )53|) into spin non-flip $2-parton n f 
and spin flip $2-parton f contributions 

^*~>-pww C ab Q 2 f , TWW , _, tww\ 

$ <Lparton,n.f. = 2 J \P A + [V ~ VWl ) 

f d 2 x 1 

x / ^ ^e;-e 1 fi\x\K 1 (fi\x\)Af(x,k) (54) 

and 

^y*-^pWW C ab Q 2 f / TWW i —\ T WW\ 
^-parton.f. = ~ J d V [<Pa ~ \V ~ Vifl ) 



x / ^ -&y ~ ' ~ ] ) ^KMx\W&Q ■ (55) 



Let us note that after the introduction of quark helicites as characteristics 
of the states (here quark and antiquark are massless, therefore the quark 
and antiquark helicities are opposite), the final formula has a form (see e.g. 
Ref. |69j) of a sum over intermediate quark helicity (A) states of the products 
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of photon wave functions, dipole factor and the specific combinations of DAs 
for each p meson helicity state, according to 

CiS. = ^J d y J d 2 %g 2 S ab J\f(x,k) 5>8^*3t (57) 

The coefficient explicitly factorised out, is related to the partonic con- 
tent of the p°— meson as a uu ~^ d state such that the meson involved below 
is understood as a one flavour quark-antiquark state. We now define the 
following combinations 



VP a +<Pi )> (58) 



1 c 



x/2 8iV ( 



v/2 8iV c 



(59) 



( ^TW_^TW ); (6Q) 



V2 8iV ( 

1 



vPa +<Pi ), (61) 



v/2 8iV ( 



c 
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Figure 5: The 2-parton contribution to the 7* — > pt impact factor in the dipole factorised 
form. 



and for the 7* [70 

®-+(y,x 



%-y 



%-y 



%-y 



-7 


X 


\x 




e+ • 
-7 


X 


\x 




^ ■ 


X 


\x 




^ ■ 


X 


1^1 



C X 

pKx{p\x\) = — j^y—pK^p^) 



[iK^nlx]) 



pK^plx}) 



7T\/2 \x\ 



71V2 k 

with, for the polarisations of the 7* and the p meson, 

/ \ 



7T\/2 \x\ 

p x * 

y—pK^p^) 



± 



=F 



V2 



1 

±i 

V / 



and 



x = X\ + ix2 ■ 



(62) 
(63) 
(64) 
(65) 



(66) 



(67) 



The obtained factorised structures (l5B"j) and (J57|) are illustrated in Fig. EJ 
We believe that the combinations of DAs appearing in these formulae are 
universal and related to the first derivative of the unknown p meson wave 
function. It would be interesting to have a deeper understanding of the 
appearance of those combinations of DAs. 

Performing the integral over x of the results (1541) and ( )55l) using 



—p—fR^p^y- 

Ztx \x\ 



S - 



(68) 
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we recover the same results in momentum space of Ref. 57J for the impact 



factors $ 



7 ~>P 

2— partem, n.f. 



and $ 



2— parton, n.f. 



2— parton, f . 

C ab Q 2 



x T n . f . 



2— parton, f . 



C ab Q 2 



xT, 



2 k 2 {k 2 + 2p 2 ) 
(^ 2 + /i 2 ) 2 ' 

dy M w - (y - y)^^) 

(fc 2 + M 2 ) 2 ' 



(69) 



with for the spin non-flip T n f and spin flip Tf tensors 

(e; -k)(k- e 7 ; 



(70) 

(71) 
(72) 



This fact can be seen as a self-consistency check of our calculation. Note 
that in the treatment in momentum space performed in Refs. 56|, [57J, the 
contribution up to twist 3 of the 2-body part involves 6 diagrams for H T>i (y), 
(a), .., (f) and 12 additional diagrams (al), (a2), ..,(f2) corresponding to the 
first order term of the Taylor expansion around the longitudinal direction p\ 
of the hard part d ^ (£ = ypi). 

2. 6. 72 — > pl impact factor in dipole picture 

We end up this section by recalling for completeness the result for the 
72 — > pi impact factor. The 72 —> Pl impact factor starts at twist 2, and the 
following result holds up to twist 3 as all contributions of twist 3 are zero. 
The only contribution involves the twist 2 correlator, 



dX 



^e- iX y{p(p)$(Xn) 7^(0)|0) = m p f p ^ (y) (e* • n)p la . (73) 



2n 



According to Eq. ([14]) . the amplitude reads 



7*~>P(0) 
2— parton 



4 

m p f p 



dX 



-~ \dy I d 2 x ± H^(y,x ± )l ^-e^ (p(j>)${Xn) 7o ^(0) |0> 



* -n) I d 2 x ± H^(y, 
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(74) 



With the same kinematics, the hard part reads 

The Fourier transform has thus the form 

M(y,x) = yyK {fi\x\)g 2 5 ab Af{x,k) 

lL \y,-> (2 7 r) v / 2 2iV c °^ui;y v_,_; 

where we identify the wave function of 

71 Q 



Kx&v) = -^K (jjl\x\) (77) 



in accordance with Ref. 7JJ (up to a normalization). The impact factor has 
thus the form 

= ^ jdy J d\ g ^U{xA) E (el 'sN c liy) ^ y) 

= ^ jdyjd\gH ab N{x,k)Y,^l\{^y) (78) 

A 

where the meson involved is again considered as a one flavour quark-ant iquark 
state, with a wave function defined as 

0oA = f^V(y). (79) 

Note that performing the integral over x we recover, as expected, the 
same result as in momentum space 

= e9 Tj: Q [ dyyyvM f ^oW)(l - e**){l - e~^) 



y/22N c J aww J 2tt 



dywx{y)- 



V22N C Q J k 2 + ^ 2 ' 
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3. The 3-parton impact factor 



3.1. Collinear factorisation beyond leading twist, in the transverse coordinate 
space 

To perform a full twist 3 computation of the ^t~^pt j m pact factor, we 
have to go beyond the WW aproximation and take into account the p— meson 
3-parton Fock state contribution. The treatment that we will follow below 





Figure 6: Factorization with Fierz identity of the 3-parton 7* — > p impact factor. 



for the 3-parton amplitude, goes globally in the same way as in the 2-parton 
case. 

Using the Fierz identity we first factorise spinor and colour indices linking 
the hard part to the fields of the non-local correlator, see Fig. 0, as 



d% dHg 

(2tt) 4 (2tt) 4 
1 f d% d 4 



tr[H a (£ u £ g )S Q (£ u £ g )} 



tT[H a (£ l7 £ g )T^]S^(£ u £ g ). (80) 



4 J (2vr) 4 (2tt) 

We use the Sudakov decomposition of the momenta of partons i, £i = atip + 
(3in + £i± and for the Fourier conjugate coordinates Zi = a Zi p + (3 Zi n + Zi± in 
the argument of the non-local correlator defining the soft part S. We factorise 
the amplitude in the momentum space, and reduce it to a convolution on the 
longitudinal fractions yi of the p meson momentum p carried by the partons. 
It reads 



x / e 

2vr 



(2tt) 4 (2tt) 4 

'-( 

2tt 



( ^lj\ 1 (e 1 -n-y 1 ) f ^9_ e i\ g (e 9 -n-y g ) I ^4 - g-«l-«l / e^ 9 '* 9 



x(p(pMz 1 ) i T f} gAi(z g )m\0) 



(81) 
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Using the Sudakov decompositions of the momenta £i and of the coordinates 
Zi we obtain the expression 

dondPx f da g d(3 g f dXidXg iXiyi _ iXgVg 



(2tt) 2 J (2tt) 2 J 2tt 27T 6 



x J da zi d(5 zi J ^ Z9 ^ g e- iQl(/3z i- Al) - iQ9( ^- A9) e- i/3ia ^- i/39a ^ 



in which after simple integrations over ctj and /?, we obtain the following 
formula in which the hard part has still a dependence on the transverse 
momenta of partons 

n-Zton = -\Jdy 1 dy g J^^tr[H a (y 1 ,y g ,£ 1 ±,i 9 ±)^} 

x J ^e~ lX ^ J d ^~ lXsVa J d 2 z 1± e- i£ ^ z ^ j d 2 z gl _e~ it ^ z ^ 
x(p(p)|^(A 1 n + z 1± ) iT p gAi(\ g n + z s± )V(0)|0) . (83) 

Let us denote by H a,r>1 (yi, yj, £i±, £j±) the trace which involves the hard scat- 
tering amplitude 

H a ' r \ yi ,y v £ i± ,£ J± ) = tT[H a ( yi , yv £ l± ,£ J± )^] , (84) 
and let us introduce its Fourier transform H T (yi,yj,Xi±,Xj±) defined as 

= J d 2 x l± d 2 x j± H a ^( Vl , Vj , x l± ,x J± ) e-*(Ax-*ix+^-* 3 -x) . ( 85 ) 

At twist 3, the 3-parton hard part contribution is given by the first term 
of the Taylor expansion around the collinear direction. Hence the 3-parton 
contribution to twist 3 is given by 

H a ' Tfi ( y i, yj ,0 ± ,0 ± ) = J d 2 x 1 ±d 2 x g _ L H a ' Td (y 1 ,y g ,Xi±,Xg_ L ) . (86) 
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Substituting the representation (155]) of the hard part in Eq. (153)) and per- 
forming the integral over we thus obtain 

®t^arton = -\ J d yi dy g J d 2 x 1± d 2 x g± H a ^ {y x , y g , x 1± , x g± ) 
^ e_iAiyi / ^e-^(p(p)|^(A 1 n)^^(V)^(0)|0).(87) 

The soft correlators appearing in Eq. (1571) are parametrised by 
d\\ dX 



(p(p)|V(Am)i 7/i ^(A s n)^(0)|0) 



2tt 2tt 

= / p (X p B(yi, y 2 ) P^ e p±Q , (88) 
dXi dX & 
2tt 2tt 



^e- iAl ^- iA «^(p(p)|^(A 1 n)i 7 57MMa(A 5 n)^(0)|0) 



/p C 3p iD(y 1 ,y 2 ) p^ e aep±pn . (89) 



Using these parametrisation in (1571) . we obtain 

= — ^ / dy x dy g j d 2 x 1± d 2 x g± 

x H e » x >*{y x , y g , x 1± ,x g± ) Cs P B (y^ Vi + Vg) 

+ H^^y^x^Xg^i^D^yx +%)] , (90) 

where we use the notations 

H a ' r ^ = H a ^a a b^ (91) 
R± a = e ailv<T e^p v rf . (92) 

It is more convenient to use the combinations S(yi,y 2 ) and M(yi,y 2 ) of 3- 
parton DAs defined in Ref. [58( as 

S(yi,y2) = ClB( yi ,y 2 ) + Ci p D( yi ,y 2 ), (93) 



:lB( yi , y2 )-ci P 

S(yi,y 2 ) and M(yi,y 2 ) satisfy the following relations 



M( yi ,y 2 ) = (lB( yi , y2 )-(£ p D( yi ,y 2 ). (94) 



S(y 2 ,yi) = -M(y u y 2 ), (95) 
M(y 2 , yi) = -S{ yi ,y 2 ), (96) 
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Figure 7: The kinematics for 3-parton contributions. 



due to the symmetry properties of the DAs B(yi,y 2 ) and Z) (2/1,2/2) under 
C— parity transformations. Using the combinations ( 19"3"| of DAs, the 
3-parton contribution f[9"Uj) reads 



= ^ ^ d 2 X 1± d 2 X g± 

.(97) 



+ M{VU f + Vg) (H^ivi, y g , x 1± , x g± ) - 1 H R ^(yi, y g , x 1± , x g± )) 



The next section is mainly devoted to the computations of the Fourier trans- 
forms H ep± ^(yi,y g ,xi ± ,x g± ) ±iH R - L ' jhB (y 1 ,y g ,Xi±,x g _L). 

3.2. Diagrams and colour structure of the 3-parton hard part 
3.2.1. Kinematics 

To get information on the dependence of the Fourier transforms H's on 
transverse coordinates, we need to extend the kinematics in order to in- 
clude the transverse degrees of freedom of the momenta of the partons. The 
kinematics we use here, as in the 2-partons kinematics Eq. ( |25|) . satisfies 
the on-shell condition for each parton, and it reduces to the kinematics of 
Ref. 57[ in the collinear limit. 

The Sudakov decomposition of the momenta of the partons are thus de- 
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fined in analogy with Eq. 

t 2 

£1 = yiPi + ii± + — P2, (98) 

Vx* 
£ 2 

£2 = V2Pi + t2± + —P2i (99) 

V2 S 
I 2 

£ 9 = y g Pi + £ g ± + —P2, (100) 

with yi + 1/2 + Vg = 1 an d £11 + ^2_i_ + = 0. The resulting momentum of 
the p meson is: 



5 

S \ Vl V2 Vg I S 



I t t t \ V 

Pp = Pl + -{- + - + -) P2=Pl + ~ L P2 (101) 



and thus 



p2 p2 f 

pi = ^ + ^ + ^. (102) 

p yi 2/2 y g 

The momenta of the incoming photon and of the t-channel gluons are fixed 
according to Eqs. (|26|) and (127|) . 

3.2.2. Classification of the diagrams in colour dipole configurations 

Below we consider QCD with SU(N C ) group, with N c finite. The fac- 
torised impact factor in momentum space reads 

fc~i*^p _ ^ m pfp f 7 j 

^3-parton,ext. ~ ^ / "2/1 

g(j^+j^ (g e fU . jt(yi) ^ ^ i±j v) + 4 H R ^( Vl , y g , £ 1± ,£ g± )) 
+ M ^^ 1 + ^ (F e ^( yi , „„ III, V) - < # *^ 75 (2/1, 4l±, V)) ■ (103) 

Note that the above expression $3_^ a rfon, ext. * s an extension of the impact 
factor, in such a way that the hard part is treated in the kinematics where the 
partons carry non zero transverse momenta. This ^Jj^arton ext mixes twist 3 
terms (which are the only one remaining in the collinear limit £±± = £ g ± = 0), 
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(cG2) (dG2) (fG2) 

Figure 8: The 12 "Abclian" (i.e. without triple gluon vertex) type contributions from 
the hard scattering amplitude attached to the 3-parton correlators for the 7* — > p impact 
factor, with momentum flux of external line, along p\ direction. 
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(atGl) 



(atG2) (btGl) (btG2) 




(ctGl) (ctG2) (dtGl) (dtG2) 




(etGl) (etG2) (ftGl) (ftG2) 

Figure 9: The 12 "non-Abelian" (with one triple gluon vertex) contributions from the hard 
scattering amplitude attached to the 3-parton correlators, for the 7* — > p impact factor, 
with momentum flux of external line, along pi direction. 



with higher twist terms introduced by the nonvanishing transverse momenta 
£i± and £ g ±. 

The panel of diagrams contributing to Eq. ( I103p is composed of 12 'Abelian' 
diagrams, shown in Fig. |8]and 12 'Non-abelian' diagrams with one gluon ver- 
tex, see Fig. [9] and 4 'Non-Abelian' diagrams with 2 triple-gluon vertices, see 
Fig. [TOj Similarly to the 2-parton hard part computation, we perform the 
integral over k, according to Eq. OH]), using the residues method applied to 
the contour C_. 



These calculations proceed in full analogy with the ones performed in (57 
For completeness of the present paper, we present in the Appendix A in more 
details the calculation of the diagrams aGl (of the Abelian types), atGl (the 
non-Abelian type with one triple gluon vertex) and gttGl (the non-Abelian 
type with 2 triple gluon vertex). The degrees of freedom which we want to 
consider now are colour dipoles, made of a colour-anticolour singlet pair in 
the fundamental representation, constructed either from a qq pair, denoted 
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(gttGl) (gttG2) (httGl) (httG2) 

Figure 10: The 4 " non- Abelian " (with two triple gluon vertices) contributions from the 
hard scattering amplitude attached to the 3-parton correlators, for the 7* — > p impact 
factor, with momentum flux of external line, along p\ direction. 




Figure 11: The {qq}, {qg} and {qg} dipoles interacting with the t— channel gluons. 

by{qq}, or from a quark and the "antiquark part" of the gluon (denoted by 
{qg}), or from an antiquark and the "quark part" of the gluon (denoted by 
{qg})- We now group, see Fig. [TTJ the diagrams involving independently 
the 3 colour dipoles constituted by the pairs {qq}, {qg} and {qg}- Due 
to the topology of the associated diagrams, the dipole {qq} is suppressed 
by 1/iVj?, the corresponding diagram being non-planar. We will justify in 
Section [3. 2. 31 that these diagrams contain only colour dipole interactions with 
the t— channel gluons in colour space. 




Figure 12: The {qg} dipole content. 
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The 6 diagrams corresponding to the interaction of the {qg} system with 
the t— channel gluons are the contribution in ^j^^f; of aGl, httGl, atGl, 
etGl, dtGl, btG2, shown in Fig.[12j The results of the diagrams associated to 
the {qg} system, cGl, gttGl, ctGl, ftGl, btGl, dtG2 are obtained from the 
diagrams of the {qg} system by exchanging the role of the quark and the anti- 
quark. The diagrams associated to the {qq} system are the contribution in 

\ Cf ~ j °^ bG2, dGl and the symmetric diagrams under exchange 
of the quark and the anti-quark, cGl, dG2, bGl, shown in Fig. d3j 
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Figure 13: The {qq} dipole content. 



3.2.3. Colour structure of the hard part 

In colour space, after factorisation of the hard part and of the soft part, 
two types of diagrams appear, the abelian diagrams corresponding to the {qq} 
system and the non-abelian diagrams corresponding to the {qg} and {qg} 
systems. Using Fierz identity we can show that in colour space both type 
of diagrams are equivalent to a colour dipole interacting with the t— channel 
gluons. For the abelian one, one gets 




(104) 



(105) 
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This second identity can be easily derived based on the relation 



Tr([t a ,t b }t c )= t -f abc , (106) 



which can be represented graphically as 



i'-= b XX c - »<U.- (107) 

thus allowing to pass from the adjoint representation to the fundamental one. 
We thus conclude from Eqs. (I104| 11051) that in colour space we expect only 
dipole colour contributions even at finite Nq- 

3.3. Fourier transforms of the 3-parton diagrams in the collinear limit 

In what follows we denote x { and £_ { respectively the transverse position 
and momentum of the parton %. In order to get connection with the dipole 
picture involving three partonic states and interacting with two t— channel 
gluons, it is important to identify the relative momentum £ { - of the partons 
forming the dipole {ij} interacting with these gluons. This will result for 
example in the appearance in the formulae (I112| 11131) of the dependence on 
the center of mass momentum of the partons forming the interacting 
dipole. As we show later, this dependence in will affect the final form 
of formulae ( II 12| II 13j) after Fourier transform with respect to £ t ■ and then 
taking the collinear limit. The 3-parton system carries two colour dipoles of 
momenta and l^. The interaction of the t— channel gluons having trans- 
verse momentum k with the dipole ij is encoded in the dependence on gluonic 
momenta having the form of the shift £^ + k. Thus in our approach only 
information about a dipole of momentum interacting with the t— channel 
gluons survives the collinear limit. 



Following Ref. [721 ] . we identify the proper definitions of £^ and in 
a way inspired by the analogy that exists in the infinite momentum frame 
between the transformations of the sub-group F of Lorentz group leaving 
invariant the hypersurface orthogonal to the direction pi and the Galilei 



tranformations in nonrelativistic 2-dimensional mechanics. In Ref. [72j, a 
dictionary is established between the generators of the Galilean transfor- 
mations of a 2-dimensional system of nonrelativistic particles of masses 
positions Xi and momenta £ i} and the generators of the sub-group F of the 
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system of particles of longitudinal fractions of momentum y iy transverse posi- 
tions Xj and transverse momenta £ i . The analog of the mass m, of the parton 
% is proportional to yiQ in our computation, where yi is respectively y±, 1/2 
and y g for the quark, the antiquark and the gluon. 



{qg} 


Center of mass G gg 


Reduced particle R gg 


momenta 


L 9B = i a + *, = -*! 


-Q9 m+y g yi 


positions 


_ V2X 2 +y g x g 
— u qg yi 


±. — ±.2 ^-g 


masses 


m Gqg = (y 2 + y g )Q 


m R _ = ^ = my 9 Q 



Table 1: Kincmatical variables of the center of mass Gq g and of the reduced particle R gg 
for the system {qg}. 

Let us focus on the diagram atGl as an example. The diagram atGl 
corresponds to the interaction of the antiquark-gluon 2-parton system {qg} 
with the t-channel gluons. According to the above discussion, we treat the 2- 
parton system {qg} as a 2-body non-relativistic system which is characterized 
by its center of mass Gq g of position XQ_ g and momentum L^ g and its reduced 
particle Rg g of relative position x from the center of mass. We show in 
Table [1] the kinematical variables associated with Gq g and Rq g for the system 
{qg}- This is illustrated in Fig. O The Fourier transform of the hard part 
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Figure 14: The dynamics of the colour dipolcs, illustrated for the diagram atGl. Here G 
is the center of mass of the 3-parton system, Gq g is the center of mass of the system qg, R' 
and R are the respective reduced particles of the 2-body system {qGg g } and {qg}- The 
blue (black) dots symbolise the color c q and anticolor c q charges, which make the spectator 
dipolc. The red (grey) dots symbolise the color Cq and anticolor Cq charges, which make 
the dipole interacting with the two t— channel gluons. 



associated to the 2-parton system {qg} reads in this kinematics 

x exp(-i((L - £ 2 - i g ) ■x 1 +£ 2 -x 2 + £ g - x g )) 
d 2 L d 2 L dH 



-2 " -9 



:H^(L,£ 2 ,£ q )S 2 (L) 



(2tt) 2 (2vr) 2 (2tt) 2 v-'^'^ 
x exp(-z(L ■ x G + L gg ■ (x G _ g - x x ) + £ gg ■ x)) 
d 2 P d 2 £ 
(27r)2(2vr)2^ 
x exp(-i(L^ g ■ (xc_ g - xj) + ig g ■ x)) (108) 

with L = £ x + £ 2 + L g an d x G — yix 1 + y<i x 2 + y g x g the momentum and the 
position of the center of mass G of the 3-parton system. We now perform 
the change of variables (£ 2 , £ ) — > (£* g , L gg ), which involves the Jacobian 
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(2/2 + Vg)/yi = 1, leading to 

m^^^hj ^ ^ ^ tel (^,^ r )exp(-i(L j& -(x 0tt -x 1 )+4 t> -x)) ■ 

(109) 

Let us denote F^y^, L^^^yj) and F^^^L^.y^yj) respectively 
the vector and axial-vector contributions of the diagram 'Diagr' related to 
the 2-parton system {ij}. They read 

Fg&y^LvMVj) = l -^l(S( yi ,y 2 ) + M( yi ,y 2 )) 

x HW^'tfay^Ly), (110) 

l-t^JLyLy'l,^) = ^ 1 -(S(y 1 ,y 2 )-M(y 1 ,y 2 )) 

x H^'^fay^Lq). (Ill) 

The sum over all the diagrams associated to the 2-parton system {ij} of 
all the F^y^L^y^yj) and F^^l^L^^y^ will be denoted by 
(lij) hijiViiVi)- Note that after using the EOM, as shown in Section |4~T1 
the {ij} 2-parton system will be identified with the dipole {ij}. 

For the diagram atGl chosen as an example, we compute the diagrams 
associated to H^i e ^(y 2 ,y g ,^ g ,L^ g ) and H^^^y,,^,^) and 
we find respectively for the explicit expressions 

FS\ y (i, 9 ,L^y 2 ,y g ) = — (S(y h y 2 ) + M( Vl , y 2 )) —- 




» + *)' 



ZymQ 

(Lfg ■ e,)(e* p ■ (4- g + k))) (L, 9 ■ ep(e 7 ■ (^ + k)) 

ZmvgQ/ (Vi + 2/2) 2?/ 2 Q 
1 

X ^T2 ^\ / r.2 ~ ^ XM2\ ' ( 112 ) 



( ^ 9 






Q 




V 2 Mi/Q 






2 
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F i«n,A(£qg>Lqg>V2>yg) = ~ ^ 7T \ ( s (vi, Vz) - M iVu Vz)) - 



2 C F 2 y ^""'^Q 



2 2y 2% Q/(y 1 + y 2 ) 



2y 2 y g Q/(yi + y 2 ) 2y 2 Q 
1 



2^7Q f 2 y ^2^7Q ^ 2 f 2^- 9 /Q 



113) 



with n\ = ymQ 2 . 

It is interesting to note that results (I112[ 11131) exhibit the appearance of 
two colour dipole dynamics, as shown in Fig. [TH 

The first dipole is made of the colour charge Cq of the antiquark and the 
conjuguated charge present in the gluon Cq. In momentum space, the colour 
Cq is at rest and the conjuguate colour Cq is carried by the reduced particle 
Rqg of the system {qg} with momentum £q g and mass rriR- g = fJ^ g /Q giving 
a dipole of size x = x 2 — x g . This first dipole interacts with the t— channel 
gluon as the momentum £ gg is shifted in + k. 

The second dipole is made of the colour charge c q of the quark and the 
anti-colour charge present in the gluon wave function c q . We can interpret 
that c q is carried by a particle at rest and c q by the reduced particle R' of 
the system constituted of the quark and the center of mass Gq g of the {qg} 
system. The particle R' has an analoguous mass 

mjf =g = (i + _^r|, („ 4) 

Q \yi y2 + y g J 2 

a momentum 

+ ».)!,- 
yi + V2 + y g 

and a transverse size x R , = x G — x 1 . This dipole does not interact with the 
t— channel gluons in this diagram and thus its contribution will disappear in 
the collinear limit. 

The kinetic energies of the two reduced particles entering the dynamics 
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of the dipoles 



E, 



c,R> 



±4g 



E r 



2/if/g ' 
2/jL/Q 



(116) 
(117) 



appear explicitely in the two last denominators of Eqs. (11121 II 1 3 j) . In the 
collinear limit up to twist 3 accuracy only information about a dipole formed 
by the reduced particle interacting with t— channel gluons survives in the 
transverse collinear space. We can then simplify the results of Eqs. (11121 
I113p . by taking the limit L^ g = 0, as illustrated Fig. [151 which gives 



F&lv(l* 9 ,Q,y2,y 9 ) 



C ab N c 1 
2 C F 2 

Vqg 



{S{y 1 ,y 2 ) + M(y 1 ,y 2 )) 



x 



^{ig} 

atGl,A\^qg 



(4 5 ,Q, y 2 ,y g ) 



C ab N n 1 



lis) 



2 C F 2 ^ Vl,v ^ ~ M ^!'^)) 



X 



/4 



yt 2(/*|, + (^ + fc) a ) ' 
We can now express Eqs. (11181 II 1 91) in terms of their Fourier transforms 



119) 




/ VW: 




-93 



-t/3 



Figure 15: The relevant momentum £^ g of the interation with t— channel gluons in the 
collinear approximation. 
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F^l v (l^0,y 2 ,y g ) = (S( yi ,y 2 ) + M(y 1 ,y 2 )) 

x ^ / Ps^K ( Ng \x\) e^+B* , (120) 
2 J {2tt) yi 

fab at i 

F { aH,A(^y2,y g ) = —^(S(y 1 ,y 2 )-M(y 1 ,y 2 )) 

/ ^^K^-M)e^--. (121) 
Finally, what remains in the collinear limit is 

fab at i 

F^{ v (0,0,m,y 9 ) = —-^-{S{y^y 2 ) + M{y^y 2 )) 

x ^r 2 / 7^— M^ 9 k\)^, (122) 

2 J (2tt) yi 

fab at 1 

^Iw)i Ko(p '' kl)e ' k -'-- (123) 

The total contribution of the diagram atGl is the sum of the above vector 
and axial contributions 

fiab at p* . p r ^72™ i,2 

The computations of contributions of all the other 3-parton diagrams go 
the same way. We first compute the diagrams associated to a dipole config- 
uration, in terms of the center of mass and the reduced particle momenta 
and masses, to obtain F^d^, L^.yi^yj). We compute then the Fourier 

transform f^(x,yi,yj) of F^(£ i: j, = 0,yi,yj) as is never shifted by 
the t— channel gluon transverse momenta k. Finally, the impact factor is ex- 
pressed in terms of f^ l ^(x_,yi,yj), where X — X j Xj IS the size of the dipole 
{ij}, i.e. the variable conjugate to the momentum The impact factor 
has thus the general form 

3— par-ton / j 3—parton \ ^ ) 

m 

= £ / d V* d Vi I i^J^^Vi)- (126) 
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The results of the diagrams in momentum space exhibit two kinds of 
structure denoted by Si(£,/j,) and S2mn{L Ha, Hb) 



Sx(£,h) 

S 2m n {L HA, Hb) 



(127) 
(128) 



where m and n are 2-dimensional euclidean indices and /i, /x^, /i# are the 
energies scales at stake. The Fourier transforms of formulae (11271) and (11281) 

are 

Si{Lh) = I ^r^K (^\x\)e^ (129) 



and 

S2mn{L HA, Hb] 

d 2 x 



1 



(27T) 

d 2 x e~ 1 -- d d 
(2tt) \i\ - n% dx m dx r 

Ornr 

Ha 



~Hb 



(2k) Ha-H 2 b\ 2 

+ Kq(h b \x\) 



K(Ha 1 


x\) 


Ha 


X 





(Mha \x\)-K {im b \x\)) 
+ K(Ha\x\) 



k o(Hb 1 


x\) 


Hb 


L 





$mn •HLm—n 



X 



Ha 



K(ha 


x\) 


Ha 


X 





~K(Ha\x\) 



'Hb 



k 'o(Hb 


x\) 


Hb 


X 





K'q{Hb \x\) 



(130) 



This expression can be simplified given that the modified Bessel function 
K U (X) satisfies the equation 

A 2 <(A) + A<(A) = (A 2 + v 2 )K u (X) . (131) 



The expression (11301) thus reads 

S2mn{L, Ha, Hb 
5, 



72 

ax e 



ii-x 



J mn r 2 



(2tt) h\ - h% 
[fi 2 A K {fi A \x\) - h% k o{Hb \x\)] 

[H A K 2 (ha \x\) - Hb K 2 {Hb 
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\x\ 



X 



(132) 



where we have used the relation 



- jK' (X) + K'^X) = K 2 (X) (133) 



implied by the standard Bessel recursion formulas [73|. Let us note that the 
Fourier transforms in Eqs. (11291 11321) lead to the appearence of two functions 
in the 3-parton impact factor, one associated with the spin non-flip transition 
ty n f and one associated to the spin flip transition ^ f\ 

* n .f. oc /j, 2 K (fi\x) (134) 
* f . oc fi 2 K 2 {fi\x) (135) 

3.4- Spin non-flip and spin flip 3-parton impact factor 

In this section, we show the results we obtain for each colour dipole 
configuration interacting with the gluons in t— channel. 

The sum of the contributions in 1^ of the diagrams (cGl), (ctGl), 
(ftGl), (httGl), (btGl) and (dtG2), associated to the scattering amplitude 
of the {qg} system on the t— channel gluons, leads to the impact factor 

CabN c f , , f d 2 x ... 

2 



■ 3-parton, {qg} 2 C F J J (2TT 



X 



M(y 1 ,y 2 ) 2 i i i \ S(y 1 ,y 2 



Vqg K o(Vqg \x\) ~ /^O^l \x\) 



qq 

V2 u yi 



Vi 

e^-e 7 e* p -xx-e 7 \ ( 'S{y u y 2 ) M(y u y 2 ) 



+ — ztM(j/i, y 2 ) [ffiKoQjt! \x\) - \i\ 2 K Q (n qg \x\)] 

ym 

+ 



x 



\x\ J \ yi 2/2 
\n\ g K 2 {n qg - n\K 2 {n x |x|)] } 



Ca " Nc I A A f d 2 %_ rp S(y u y 2 ) 2 

+ -ittt / dyidy 2 / -7^T n . L = ^0(^1 \x\) , (136) 

2 C F J J (2it) yi 

2 

with H2 = y2l/2Q 2 - Note that ^ = y 2 y 2 Q is associated to the analoguous 
reduced mass of the 2-body system constituted by the antiquark and the 
center of mass of the quark and the gluon. We show in the Table HJ the 
kinematical variables associated to the center of mass G and the reduced 
particle R of the system {qg} that we use to obtain, after simplifications 
described previously, the result 
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{qg} 


Center of mass G 


Reduced particle R 


momenta 


L-qg — L\ + Lg — ~i.2 


11-, p 1/ p iii-. p 11 p 

£ tfl i 9 Vg Li yl Lg Vg £1 

-19 ~ yg+yi y2 


positions 


Jb r~i — 

— (_r y 2 


Jb — y 1 


masses 


m c = yiQ 


mR = & = yiygQ 



Tabic 2: Kinematical variables of the center of mass G and of the reduced particle R of 
the system {qg} 



The result for the {qg} dipole is straightforward by exchanging the role 
of the quark and the antiquark i.e. exchanging y\ and i n ( H36p . Adding 
the results for the {qg} and for the {qg} dipoles and using the symmetry 
properties of S(yi,y 2 ) and M(yi,y 2 ) given by Eqs. ( 195]) and f )96|) . the spin 



non-flip part ^^ {qg}+m 



and the spin flip part ^an m ,{i9}+m read 



( | ) 7*^p,n.f. 

3-parton,{qg}+{qg} 



-C ab N, 



C 



2 Ct 



x 



+ 



X 



dy 1 dy 2 



d 2 x 



[t4 K 0(Vl \X\) + fifgKoblgg \X\)] 



£(2/1,2/2) 
yi 

M(y 1 ,y 2 ) r 2 , . h 2 . . 

/i 2 A (At2 \X\) + K (Hgg \X 

2/2 

[H 2 Kq(h 2 \x\) - n qg K (^g \x 

- [fijKoifii \x\) - fi 2 qg K (fi qg \x 
d 2 x e* ■ e 7 



2/13/2/ yi 

2/22/1 \ M(yt,y 2 



ym 

C ab N, 



c 



2 C F 



2/2 

dyidy 2 



(27T) 



^(2/1, 2/2) 1 h M(j/i,y 2 ) 1 1 \ 

ViK (ni \x\) /i 2 A o(/^2 k|) 



2/i 



2/2 



(137) 
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and 



3-parton,{qg}+{qg} 



C ab N, 



C 



2 C F 
S(yi,y 2 ) M{y u y 2 ) 



dyidy 2 



d 2 x 



e p -e 7 e p -xx-e^ 



yi 2/2 
+f4g K ^9 k\) ~ A K 2^2 \x\j\ . 



(27r)""^'^ /v 2 
[vlgK^qg k\) - f4K 2 {fJ,i \x\) 



\x\ 



(138) 



The spin non-flip and spin flip impact factors associated to the scattering 
amplitude of the dipole {qq} on the t— channel gluons are given by the con- 
tributions of type 2) from the diagrams (aGl), (cGl), (bGl), (dG2), 
(bG2), (dGl). The results read 



<| ) 7*->P > n.f. 
3— parton, {qq} 



c 



ab 



N, 



c 



C f 



dyidy 2 



d 2 x 

W) 



N(x,k) 



S(yi,y 2 ) i in M(y 1} y 2 ) 2ir , , h 

fi^obii \x\) V 2 K {fi 2 \x\) 

yi 2/2 
S(yi,y 2 ) 



y 9 



— {nlK Q (fl 2 \X_\) - fi 2 qq -Ko(flqq \X\)) 

1/2 



+ M ^ VUV ^ [(fllK {fi 2 \X\) - fl 2 q gK Q (figg \X_\)) 



V 9 



+ 2/i ^ lKo ^ 1 '-'^ _ V 2 m K ^m I^W 



c 



ab 



- 2 



d 2 £ e* • e 



±. — p —7 



(2 



7T 



S(yi,y 2 ) 2 , \ w M (2/i>2/2) 2^/ in 
AtiA (A*i |#|) ^2^0(^2 



and 



7* -»•/£>, f. 



2/i 



C ab /iV, 



2/2 



(139) 



3— parton, {qq} 2 

5(2/1,2/2) M(y h y 2 ) 



d 2 x 

W) 



N{x,k) 



^Lp ' ^L'y ^Lp X_ X_ 



\x\ 



X 



2/i 



2/2 



— {^qqK^qq ~ ^K 2 {^ \x\)) 

Vg 



■— {Aq K <l{llqq |z|) ~ l4 K 2^2 \x\)) 

Vg 



(140) 
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We show in Table [3]the kinematical variables associated to the system {qq}. 
The total 3-parton results for the spin non-flip amplitude is thus 



{<?<?} 


Center of mass G 


Reduced particle R 


momenta 


L qq = u + t* = -4 


£ _ 52^1-2/1^2 _ mli-yil 2 

-qq yi+y2 y g 


positions 


_ V1X.1+ V2X 2 
— u Vg 


Jb 1 Jb 2 


masses 


m G = y g Q 


rrll _ >4q — yimQ 
~ Q ~ y g 



Table 3: Kinematical variables of the center of mass G and of the reduced particle R of 
the system {qq} 
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(jj7*->P>n-f- 
3— par ton 



c 



ab 



dy 1 dy 2 



d 2 x 

W) 



+ 



Nc 
C F 

+ 



^(2/1, 2/2) 2 j, f I i\ M(y 1 ,y 2 ) 2 

Ml^Wl kl) ^2 K 0{^2 \X\) 

5(2/1,2/2) 2 



1J2 



Vi 



fig gKo(Hqg \X\) 



M( yi ,y 2 ) 



>J2 



H 2 a K Q (n qg \x\) 



2/22/1 
2/22/1 



X 



5(2/1 .2/2) 



2/1 

M{y u y 2 ) 



[filK (fi 2 \x\) - fI 2 gg K (fiq g \x 



+ 



TV, 



- 2 



2/2 

5(2/1,2/2) 



[/x?^o(Mi N) - n 2 qg K Q {n qg \x\)] 



y 9 



\x 



+ 272 ^ K °^ 2 '-'^ ~ V 2 qq - K ^qq \X 
([lI 2 2 K (fi2 \X\) ~ fI 2 g5 K (fX q q \x\)] 



M( yi ,y 2 ) 
Vg 



+ 2/i ^ K °^ 1 '-'^ ~~ ^lq K °^iq k 



+ 



ab 



^2/1^2/2 



5(2/1,2/2) M(y 1 ,y 2 ) 



while the spin flip 3-parton impact factor is 



<ib 



3— partem ^ 



dyidy 2 



d 2 x 

W) 



M{x,k) 



U2 



(141) 



(et • x) (x • e. 



5(2/1,2/2) M(i/ b !/2)\ JiV C , 2 2 

— J { cv ™ 99 ~~ ^i^ 2 ^ 1 Is) 



+ 



Vi 

Nc 
C F 



+^ 2 qgK2^qg \x\) ~ ^2^2 \x 



- 2 



— {v 2 q q K 2(Vqq \x\) - HlK 2 {V\ \X\)) 



Vx 
Vg 



(lI 2 qq K 2 (lIq q |x|) - [4K 2 ([1 2 \x\)) 



(142) 



In the formula (I14ip . the last line is not proportional to the dipole factor 
N(x, k). In the following part, we will show that putting together the 2- 
parton result (beyond WW approximation) and the 3-parton result, all parts 
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of the impact factor which do not have the dipole form cancel each other 
using the QCD EOM. This will extend to the full twist 3 result, the reasoning 
leading in the WW approximation from Eq. (150]) to Eq. (153]) . 



4. The twist 3 7^, — )• pr impact factor in the dipole picture 

4-1. The dipole picture arising from the equation of motion of QCD 

The relation between DAs due to the QCD EOM, at full twist 3 (it extends 
the Eq. fT49|) beyond the WW approximation) reads 

dy 



yy 



(2^3(1/) + (y - y)<Pi(y) + <fA(y)) 



+ dy ± dy 2 



S{yi,y 2 ) M(y 1 ,y 2 ) 



IJi 



y-2 



0, (143) 



with (p(y) = <fi ww {y) + (p 9en {y) being the complete DAs, i.e. which include 
both the WW and the genuine twist 3 contributions. The 2-parton impact 
factor fl50l) before using the relations due to QCD EOM reads, 



C ab Q 2 , 



<f,7 ^P _ 

2— par ton 



d 2 x 

Z7T 



x ( [(v - y)<fl{y) - <pl(y)} §* P ■ + ^ T A {y)^e* p ■ e 



p -7 



% [(?yy vsiy) + {y- y)<pT + Va(s/))] 
yy 



(144) 



Collecting all terms arising from Eqs. fll4ip and ( I144p which are not propor- 
tional to the dipole factor, we observe that they cancel 



C 



ab 



dy 

yy 



+ dyi dy 2 



yi 2/2 

due to the relation (j!43j) . The final 2-parton impact factor is thus 
QabQi r r A2. 



0,(145) 



$7 -+P 
2— parton 



f d 2 x 

dy j 2^ { [iv - v)<pT(v) - ^(y)] 



x (g;-g)(yg 7 ) +y r (y) ■ |g 



(146) 
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and it can be decomposed into the spin non-flip and the spin-flip parts 

_C ab Q 2 
2 

d 2 x 1 



2— parton, n.f . 



dy {<pa + (y- y)v>T) 



X 



2?r 2 



e* P ■ e 7 /i|x|^i(yu|x|)A/'(x,fc) 



(147) 



and 



2— parton, f . 



C ab Q 2 



d 2 x fl 



2ix \2~ p '-^ 



dy (va -(y- y)<pT) 



\x\ 



n\x\Kih-i\x\)N(x,k)- (148) 



The results f lTM [T471 IT48D are the extension of the formulae fl53l 1541 1551) to 
the full solution of the QCD EOM for the DAs, including the genuine twist 
3 solutions. 

The 3-parton spin non-flip result after using the relation (j!43p is thus 



( j ) 7*->p,n.f. 
3— parton 



X 



C 



ah 



dy 1 dy 2 



d 2 x 



N{x,k) 



Nc 
C F 

+ 



S(y h y 2 ) 2l s I i i\ M (yuy2) 2ir , i in 

= VlKo(Hl \X\) ^ 2 K {fI 2 \x\) 

yi 

~S(yi,y 2 ) 



<J2 



M(y u y 2 ) 



!J2 



mm 
myi 



S(yi,y 2 ) 
m 

M( yi ,y 2 ) 



[n 2 2 K Q (fi 2 \x\) - fi 2 qg K (fig g \x\)] 



N, 



o 



2/2 

S{y\,y 2 ) 



[li\K Q {Hl \X\) - fl 2 qg K {figg 



\x 



Vg 



{[l4 K 0{Hl \X\) - nlg-KoiUqq \X\)] 



m 
m 

M( yi ,y 2 ) 



+ — [HIK0{H2 \X\) - H 2 q K (fI q g \X 

y2 



Vg 



([nlK (n 2 \x\) - n 2 q gK (fi q q 



\x 



m 
yi 



[t4 K 0(Vl \X\) - nlgKoiflqq \X\)] 



, (149) 
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and it can be rewritten in a more compact way, using the symmetry properties 
of the DAs under exchange of y% and y 2 , as 



( j ) 7*->p,n.f. 
3— partem 



1 



-c 



ab 



dyidy 2 



^(x,k)^S( yi ,y 2 ) 



N, 



< - \2f4KoQi! \x\) + ^ [fj^KoQigg \x\) 



+ 



1 (N, 



y 9 \c 



> c 



[fi 2 2 K {fi 2 \x\) - nlgKoifigg \x 



2) [[^K (^\x\) - ^K ( Nq \x\)} 

+ — [t4 K o(V2 \x\) ~ H 2 qq K (Hqq \x\j\ 



(150) 



4-2. Equivalence with the results performed in momentum space in the light- 
cone collinear factorisation scheme 

The integration of the spin non-flip result over x is straightforward by 
using the relation 



d 2 x 



2k 2 



(151) 



The integrated result thus gives 



= ~C ab I dy 1 dy 2 e;-e 1 S(y 1 ,y 2 ) 



x 



2 - 



c 



a 



Vi 



V2 2/2 



C F J y g a + yiy 2 \a + y x y\ a + 2/22/2 
Nc oi a 2 a 



C F ay 1 + y 2 y g a + y 2 y 2 y±a + 2/12/1 



(152) 



with a = k 2 /Q 2 . The result <$M) is, as expected, identical to the one 
obtained in Ref. [57| using the light-cone collinear factorisation. 

For the spin flip result, using the symmetry of the amplitude under the 
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exchange of y\ and y 2 , we get 

dyidy 2 



x 



( J ) 7*->P, f. _ S~1Q 

3— partem 



S(yi,y 2 )\ f N c 



d 2 x 

W) 



(e!-x)(x-ej 



+ 



2/i 

Nc 
C F 



[l4g K ^i9 Isl) ~ l4 K 2(Vl \x\) 



+fi 2 g g K 2 (^ g \x\) - i? 2 K 2 {[i 2 |z|)] 

— (^qq K ^qq \x\) - l4 K 2^1 \x\)) 

Vg 

+ — {^qq^^qq ~ /A K 2^2 \x\)) 

Vg 



. (153) 



We now integrate over x = (|x| cos(9), \x\ sin(9)), with A; = cos(4>), \k\ sin(4>)) . 
Using the fact that only the spin flip contributions are non zero, and based 
on the following identities 



[{e-y -x] [x-e+] 



\x\ 



\x\ e 



V2 



[{e+Y ■ x] [x ■ e~ 



\x\ 



\x\ e 



-iO 



V2 



\x\ e 



\x\ e 



10 



\X\ 



V2 



J26 



-i2Q 



(154) 



\x\ 



(155) 



resulting from the definition of the polarisations in Eq. f )66|) . we obtain 



d 2 x 

(27T) 



Af(x, k) n 2 K 2 (fi \x\) 



de 



((eZy-x)(x-e±) 



x 



— / d\XK 2 (X) I — 2 ( 1 - cos 



kX 



cos 



J dXXK 2{ X) J d ±2 
l -e^2 J dXXK 2 {X)J 2 {^\ 



1 — cos 



kX 

fJL 



cos(#) 



„±i28 



((ejy-k) (k-e 



2k 2 
k 2 + H 2 



(156) 
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with A = \x \x\. Hence the spin flip impact factor integrated over x reads 



-2C ab T f j d yi dy 2 S( yi ,y 2 ) 



N, 



o 



ay* 



C F 

+ U~ 2 



Vi 



+ 



{)■! 



« + ym)(y20i + yiy g ) (a + y 2 y 2 ){y x a + y 2 y g ) 
oiy% ( yi y~2 



y g a + yiV2 V« + 2/12/1 « + 2/22/2 



(157) 



which is the same result as the one obtained in Ref. 57 



4-3. Complete twist 3 result of the 7^ — > px impact factor 

Combining all the 2-parton and 3-parton results for the spin non-flip and 
spin flip impact factors $^ f _5>p , $^ of the 7^ — > px transition, we finally 
obtain 



x 



^7*->/9,n.f. 

Q 2 



-c 



ab 



d 2 x 



M{x,k) 



(2tt)"^^ 2 
dy {vat{v) + (2/ - V)<Pit(v)) mN-^iG"|z|) 



12 S( " yi,y ^ [2fi 2 1 K (fi 1 \x\) + ^ [nl g K ([iq g \x\) 



Nc , 

2/1 Cf 



2/2 2/1 
2/2 2/1 



[nlK Q {n 2 \x\) - i^gKoifiqg \x\)] 



+ — (7^ - 2 J [[pilKo{fii \x\) - filgK (pi qq \x\)] 

y g j 



+ — [/^o(^2 - /J, 2 qq -K (ll q q \x\)] 

2/2 
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and 



ab 



d 2 x 



N{x,k) 



eye, 



(e*-x) {x-eS 



x 



2 



(27r)'-^'^V 2 m 
dy {<Pat(v) ~(y- y)<Prr(v)) v\x\Ki{p\x\) 



dy\dy 2 



s{yi,V2) 



Vi 



— - [n 2 g K 2 {fi qg \x\) - ffiKsQn \x\) 



Nc 
C F 



-^lg K -i{^qg \x\) - ll\K 2 (lH \x\)] 



— {^qq K 2^qq \x\) ~ /^l^Ml \x\)) 



+ — {^qq K 2^qq \x\) ~ (%K 2 (p2 \x\)) 



(159) 



The Eqs. ( 1158P and fl!59[> are the main results of the present paper. We 
have achieved the goal of presenting the 7£ T — >■ impact factor, in the 
forward limit, in a factorised form involving explicitly the coupling of two 
t— channel gluons with a dipole of transverse size x produced during the 
transition of the 7* to the p— meson. The colour structure, after factorising 
a global normalization 1/(N 2 — 1), has the form 



6 ab 



(2N C )(N 2 -1) 



[-^planar N c -\- v4 non _ planar] 



(160) 



where A p i anar and A non _ planar are both N c independent. 

This form of the impact factor will permit to include effects of saturation 
in the proton for exclusive p meson electroproduction both for longitudinally 
and transversally polarized meson 74|. The <^ s defined by Eqs. (1561 [571) 
could be extended beyond WW approximation, using Eqs. (11581 11591) . For 
this we would need to disentangle the 3-body wave-function of the trans- 
versely polarised photon 7^ , which, to the best of our knowledge, is unfor- 
tunately unknown. 



5. Conclusion 

In this paper we have shown that the ^* LT — > pi,T transitions, which 
involve contributions up to twist 3, preserve the usual dipole form of the 
scattering amplitude at finite N c . Such a dipole picture of the scattering 
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process has been used successful to describe phenomenologically the effect of 
gluonic saturation inside the target at asymptotical energies for inclusive as 
well as for exclusive processes involving twist 2 contributions treated in the 
collinear approximation. 

At first sight, twist 3 contributions either of kinematical origin or due to 
genuine 3-parton correlators have a form which does not exhibit any dipole 
structure. It is only after the use of QCD equations of motions that all terms 
can be reorganized in a way consistent with the dipole picture. In particular, 
our result shows that no quadrupole contribution is involved in this process 
at the twist 3 order. We expect that this remains valid also beyond the 
twist 3 approximation. The two contributions which we have obtained in 
Eqs. ( I158p and (j!59p have the form (I160p of a planar (scaling like N%) and a 
non-planar topology (scaling like 1). 

The obtained formula are a starting point for further investigation of the 
effect of gluonic saturation for exclusive processes beyond leading twist. 
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Appendix 

The contributions of 3-parton correlators are of two types: the first one 
being of " Abelian" type (without triple gluon vertex, see Fig. [S]) and the sec- 
ond involving non- Abelian coupling with one triple gluon vertex (see Fig. [9]) 
or two (see Fig. [TO]) . Let us first consider the "Abelian" diagrams. They 
involve two kinds of Casimir invariants: 

1 8 ab () ah 

— Tr (t c t a t b t c )=C F — = C a — : (aGl), (cGl), (eGl), (fGl) (161) 
1 / /V\ h ab h ab 

w Tr{^^)=\C F -^ — =C b — : (bGl), (dGl), (aG2), (cG2), 

(bG2), (dG2), (eG2), (fG2), (162) 

where the 1/N C comes from the Fierz coefficient when factorising the quark- 
antiquark state in colour space. To illustrate the method, we consider the 
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Figure 16: The detailed structure of the diagram (aGl). 



particuliar diagram (aGl) of Fig. |8l illustrated in Fig. [16] The vector con- 
tribution reads 



1 

raft 



12 Slf 
®aGi = -e q ^-(^9 2 f P m p jj^— dy 1 dy 2 B(y 1 ,y 2 ) (163) 

o 

dn Trfa (h - 4) f T (A +t g -4) h (fa + h) h M 



2tt [(4 - q) 2 + irj] [{h + £ g -q) 2 + irf\ [(k 2 + £2? + irj\ 



and it equals 



i 

-,2 tab 



®aci = ^y- f P m P jjf J d Vi d V2 B(yi, y 2 ) 



x { e* p • e 7 (y 2 £\ -y 1 £ 2 -£ 1 + yi y 2 Q 2 ) + y x ({£, ■ e 7 ) (£ 2 ■ e* p ) + {£, ■ e* p ) (£ 2 ■ e 7 )) } 

x fi64) 

(£% + nl) (y g V 9 ll + 2 yi y g £ x ■ £ g + ym £ 2 g + yi y 2 y g Q 2 ) ' 

The expression (j!64p is used in the description of two dipole configurations: 
{qg} and {qq} types. The changes of variables defined either in Table [U for 
the {qg} system or in the Table [3] for the {qq} system leads respectively to 
explicit forms for F^ v (£^ g , L^ g , y g , y 2 ) and F$$ v (£^, L^ q , y u y 2 ), defined by 

Eq. (rrroji . 
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The corresponding axial contribution reads 

1 

i 2 5 ab 1 f 

®aGi = - e q^-(^)9 2 f P m p ^— d yi dy 2 e a e , Tpn D(y u y 2 ) (165) 

o 

dn Trfa (fa - 4) 7a (fa + fa - 4) h (fa + fa) h M 



2tt [(£ x - q) 2 + in\[{h + £ 9 - q) 2 + irf\[(k 2 + h) 2 + iv] 



and equals 



i 

,2 sab 





x {e; ■ e 7 (y 2 £ 2 + y 1 £ 2 -£ 1 + Vl y 2 Q 2 ) - Vl [(£, • e 7 ) • e* p ) + (£, ■ e* p ) (£ 2 ■ e 7 )]} 



(166) 

(ll + tf)(y 9 y 9 ll + iviVgk ■ 4 + 2/12/1^ + ymy g Q 2 ) ' 

Similarly to the vector contribution (I164p . the expression (I166p is used again 
in the description of two dipole configurations: {qg} and {qq} types. Again, 
the changes of variables defined either in Table [U for the {qg} system or 
in the Table |3] for the {qq} system leads respectively to explicit forms for 
F^M^L^m) and F^ A (£_ qq , ^, Vl , y 2 ) , defined by Eq. (TUID. 

Consider now the " non-Abelian " diagrams of Fig. [91 involving a single 
triple gluon vertex. They involve two kinds of colour structure: 

9 /V 1 fr ab 

- % )Tr{t c t h t d )r d = ^—— : (atGl), (dtGl), (etGl), 



N 2 -V ' v ' J 2 C F 2 N c 

(btG2), (ctG2), (ftG2), (167) 

9 N 1 b ab 

' -i) Tr{t c t d t b ) f cad = — -— : (ctGl), (btGl), (ftGl), 



N 2 -V ' v ' J 2 C F 2N C 

(atG2), (dtG2), (etG2) , (168) 

where the 2/(iV? — 1) comes from the Fierz decomposition when factorising 
the quark-ant iquark gluon state in colour space. Let us consider the diagram 
(atGl) of Fig. ij illustrated in Fig. HTJ We denote as 

h u n p 4- k p n u 

d»<>(k) = g v <> - , (169) 
k ■ n 
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Figure 17: The detailed structure of the "non-Abelian" (with one triple gluon vertex) 
diagram (atGl). 



the numerator of the gluon propagator in axial gauge, and 

fc 2, h) = (h - fc 2 ) w g^ 2 + ■■■ (170) 

the momentum part of the 3-gluon vertex, where ki are incoming, labeled 
in the counter-clockwise direction. The contribution of the diagram (atGl) 
then reads, for the vector DA, 

v 12 (-i)N c 2 5* 1 f 

o 

■^ e ^.(171) 



d^ih-QV^i-h+ig, k u -l g ) 



[{h - qf + ir]][{ki - t g y + i v ][(k 2 + t 2 f + i V ] 

Note that for this diagram, as well as for all "non-Abelian" diagrams, one 
can easily check that only the g vp part of ( 11 69[) contributes. 

Closing the k contour above or below gives for the vector DA part of the 
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diagram (atGl) the result 

,2 tab 



p n 2 h ab N C 

Cgi = ^frn p f p ^ w --^jdy 1 dy 2 B(y 1 ,y 2 ) (172) 



o 



x {yi 4 • e 7 (y g £ 2 ■e;-2y 2 £ g -e* p + (ft + y 2 )k- e* p ) 

+ Vg [e* p ■ e 7 (y 2 £j -yxi x - {l 2 +k)+Vi V2 Q 2 ) + J/i & • §* {£ 2 + k) ■ ej } 

yi 



X- 



(ll + vi)(y2y 9 l 2 +yiy 9 l 2 +yiM g +yiyik +2yi{y a l 2 - ml g ) -k + yiygmQ 2 ) 

Similarly, the contribution of the diagram (atGl) reads, for the axial DA, 

, % 2 (-i)N e 2 , 5 ab 1 f , , 

atGi = ~e q 1--2CT9 m pfpJW c Ys j dVl dV2 D{Vl > m) 



x j^Tr[^(t 1 -i)^^ 2 + t 2 )^ ll5 ] 

d^(h - e g ) v P x a (-h + e g , h, -Q 

[(i 1 - q )2 + %r] ][(k 1 -£ g ^ + t V }[(k 2 + £ 2 y + i r ] } P2 1 1 

and closing the k contour above or below gives 

1 

p a 2 <) ab N f 

= - J frn p f p —-^-Jdy 1 d y2 D( yi , y2 ) (174) 



x {yik • e* (y g £ 2 ■ e 7 - 2y 2 £ g ■ e 7 + (f/i + y 2 )k- ej 
- (ep ■ e 7 ) [2/1^ ■ (y g £ 2 - 2y 2 £ g ) + ^(ft + y 2 ) k ■ £ x + y 5 y 2 ^ + yi^Z/gQ 2 ] 
[di ■ e* p ) [e 7 ■ (y g £ 2 - 2y 2 £_ g + (ft + + • e 7 )[e; • (| 2 + fc)]] } 

2/1 



x 



(£1 + (*i)(V2y g & 1 +ViVg£2+ ymlg + ymk + 2y 1 {y g £ 2 - y 2 £ g ) ■ k + yi y g y 2 Q 2 ) 
The expressions (jl72j) . (1174|) . supplemented by the changes of the variables 

^99' —99' 



defined in Table (U lead to the explicit forms for F^q\ v (£^ ,y 2 ,y g ) and 



FatGiAkgi^y^yg) ' S iven respectively by Eqs. flTTjP and (JIT3). 

We consider now the "non-Abelian" diagrams of Fig. [TUl involving two 
triple gluon vertices. They all involve the colour structure 

9 N h ab 

Tr[t c td]r a f e db= cO_ (175) 



iV c 2 -l 1 IJ J C F 2N t 
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For illustration, let us consider the diagram (gttGl) of Fig. [TOj illustrated in 
Fig. [13 It reads for the vector DA, 



i 



v 12 (-i)N e 2 5 ab 1 f 



x f ^Tr[^ (A - i) 7 „£] d»»(-q + l x + £ 2 ) 



x 



V pXa (q -it- £ 2 , h, -k 2 - £ g ) d a ?(k 2 + £ 



1)1 



[{h -q) 2 + ir l }{(-q + £ l + £ 2 ) 2 + i V ] [(fc 2 + £ g ) 2 + irj] 
x Vp T s(k 2 + £ g ,-k 2 ,-£ g )p$p T 2 e* T 5 . (176) 

It equals, when closing the k contour below on the single pole coming from 
the third propagator, 



x {2y 1 (£ 1 • e 7 ) {£_, + £_ 2 ) ■ e* p + y g (e* ■ e 7 ) + Vl Q 2 ) } 

ym 



X 



ill + tf) (ym ll + ym ll + 2ym k-£ 2 + ymy g Q" z 



The expression (11771) is used in the description of the dipole configuration 
{qg} ■ The changes of variables defined in the Table [2] lead to the explicit 
form for F^ iy (£ qg , L qg , y u y g ) , defined by Eq. (THUD. 



55 



The axial DA contribution from the diagram (gttGl) reads 

*a i2{-i)N c 2 6 ab 1 f 

* 9 ugi = - e «r s -^-9 m p f p ^ w --jdy 1 dy 2 D(y 1 ,y 2 ) 

o 

x I ^Tr[^ (A - i) 7 „£ 75 ] <r'{-q + h + £ 2 ) 

V P x a (q -h- i 2 , h, -k 2 - i g ) d^(k 2 + £ g ) 

x [{h - q) 2 + H [(-<? + h + £ 2 ) 2 + H [(h + £ g ) 2 + iv] 

x Vprrfa + £ g , -k 2 , -£ g ) V \v\ e% pn . (178) 

It equals, when closing the k contour below on the single pole coming from 
the third propagator, 

1 

p n 2 t) ab N 1 r 

Kugi = m p f P J dyi dy 2 D(y u y 2 ) (179) 



x {2 Vl (£, ■ ep & +£ 2 )-e ry -(e* p - e 7 ) [2y x £_,-£, + ( Vl + y 2 )£ 2 x + yi y 9 Q 2 ] } 



X 



(ll + t4)(ym ll + vm ll + 2ym k-l 2 + ymy g Q 2 ) ' 

Similarly like in the case of the vector contribution, the expression (11791) is 
used in the description of the dipole configuration {qg}. The changes of vari- 
ables defined in the Table|2]lead to the explicit form for F^q^^, L ?5 , y x , y g ) , 
defined by Eq. ffTTTf . 
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